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ABSTRACT 


This paper critically reviews available 
prediction methods and provides an assessment of 
their strengths and weaknesses. The methods are 
applied to selected problems which represent the 
major aero/propulsion interactions for short 
takeoff and vertical landing (STOVL) aircraft of 
current interest. The first two problems deal 
with aerodynamic performance effects during 
hover: a) out-of-ground effect, and 

b) in-ground effect. The first problem can be 
evaluated for some multijet cases; however, the 
second problem is very difficult to evaluate for 
multi Jets. The ground-environment effects due 
to wall jets and fountain flows directly affect 
hover performance. In a related problem c) hot- 
gas ingestion affects the engine operation. 

Both of these problems as well as jet noise 
affect the ability of people to work near the 
aircraft and the ability of the aircraft to 
operate near the ground. Additional problems 
are d) the power -augmented lift due to jet-flap 
effects (both in- and out-of-ground effects) and 
e) the direct jet-lift effects during short 
takeoff and landing (STOL) operations. The 
final problem is f) the aerodynamic/propulsion 
interactions in transition between hover and 
wing-borne flight. Areas where modern computa- 
tional fluid dynamics (CFD) methods can provide 
improvements to the current computational cap- 
abilities are identified. 

FOR OVER THIRTY YEARS, vertical or short takeoff 
and landing (VSTOL) aircraft have been developed 
primarily from experimental investigations 
because there are many unique aerodynamic/ 
propulsion-induced effects which occur with 
these vehicles. These problems are especially 
significant at low speeds from hover to wing- 


borne flight. Proper design of VSTOL aircraft 
requires an accounting of these complex phenom- 
ena. Unfortunately, the flows involved are not 
amenable to purely theoretical predictions at 
this time. Improved prediction methods are 
needed to guide experimental programs and to 
reduce the testing requirements in the future, 
as well as to provide an understanding of the 
flow physics. One motivation of this paper is 
to identify areas where modern computational 
fluid dynamics (CFD) methods can provide 
improvements to the current computational capa- 
bilities. The ultimate goal would be to develop 
the ability to completely analyze the aerody- 
namics throughout the entire flight envelope of 
an arbitrary VSTOL aircraft. 

Traditionally, prediction techniques appli- 
cable to VSTOL aircraft (1-8)* consist of 
potential-flow methods with empirically derived 
corrections and/or adjustments to account for 
viscous effects which significantly affect 
aerodynamic/propulsion interactions. Some of 
these methods are simple to use but provide 
results with only limited applicability. Other 
procedures build on potential-flow-panel 
methods. The input data tend to be complex and 
difficult to prepare, and often the results have 
limited usefulness. 

In recent years there has been rapid pro- 
gress in the development of CFD. This progress 
is due in part to improvements in digital com- 
puters (9) which provide rapid computation 
speeds (up to hundreds of millions of floating- 
point operations per second), and larger memo- 
ries (up to hundreds of millions of words). 

This growth was clearly described in Jamesons 
survey (10) of the evolution of computational 


•Numbers in parentheses designate 
references at end of paper. 
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methods for conventional aerodynamics. Numeri- 
cal approximation methods were examined for a 
hierarchy of models in ascending order of com- 
plexity, ranging from the linearized potential- 
flow equations to the Reynolds-averaged Navier- 
Stokes equations. In a study (11) of the influ- 
ence of CFD on experimental aerospace facili- 
ties, four main stages (I-IV) of approximation 
to the full Navier-Stokes equations were identi- 
fied. For completeness in relation to short 
takeoff and vertical landing (STOVL) aero/ 
propulsion interactions, the present authors 
have added a Stage 0 to represent the empiri- 
cally based methods, and Stage lb to represent 
linearized inviscid methods which include an 
approximate model for jets: 

Stage Approximation 

0 Empirical methods 

1 Linearized inviscid 

Ta Linearized inviscid plus boundary 

layer (BL) 

lb Linearized inviscid plus modeled Jets 

II Nonlinear inviscid 

I la Nonlinear inviscid plus interacting BL 

III Reynolds-averaged Navier-Stokes 

IV Full Navier-Stokes (large-eddy 
simulations with small-scale 
turbulence modeling) 

Stage 0 has been the mainstay method for extra- 
polating experimental data beyond the original 
wind tunnel test conditions up to the present 
time for VSTOL aero/propulsion interactions. 
Stage I is currently used for engineering anayl- 
sis of conventional airplanes. The Navy V/STOL 
Aerodynamics, Stability and Control Manual (4,5) 
utilizes Stages 0, I, and lb methods, and repre- 
sents the current state of the art in V/STOL 
computational methods. For conventional air- 
planes without propulsion- induced effects, 
limited to moderate use is also being made of 
stage II solutions; Stage III is expected to be 
widely used in the next 10 years. The principal 
items pacing introduction of Stage III into the 
aerodynamic design process are: a) development 

of improved turbulence models; b) geometry grid 
generation difficulties and limitations; 
c) development of more efficient and reliable 
numerical algorithms; and d) development of more 
powerful scientific computers. 

At this time some initial efforts have been 
made to use stage II and III approximations for 


VSTOL. One of the first stage III analyses was 
by Baker (12). This calculation of the flow 
field caused by a jet in a crossflow used a 
finite-element, pressure-interaction, parabolic 
approximation to the Reynolds-averaged full 
Navier-Stokes equations. The other papers 
(13-16) in the present session will present 
current examples of several Stage II and Stage 
III approximations applied to either ST0L- or 
STOVL-aircraft-related flow fields. 

This paper critically reviews available 
stage 0 and stage I prediction methods, and 
attempts to provide an assessment of their 
strengths and weaknesses. Existing methods will 
be discussed based on six selected problems 
which represent the major aero/propulsion inter- 
actions for STOVL aircraft of current inter- 
est. This paper presents examples using avail- 
able methods. The first two problems deal with 
aerodynamic performance effects during hover: 
a) out-of-ground effect and b) in-ground 
effect. Both problems are reasonably well 
understood for the case of a single jet. The 
first problem can be evaluated for some multijet 
cases; however, the second problem is very 
difficult to evaluate for multi jets. The next 
three problems are related to operations near 
the ground: c) the ground-environment effects 

due to wall Jets and due to surface breakdown, 
d) hot-gas ingestion effects on the engine 
operation, and e) short takeoff and landing 
(STOL) ground effects. These problems affect 
both the ability of people to work near the 
aircraft and the ability of the aircraft to 
operate near the ground. The final problem 
f) is the aerodynamic/propulsion interactions in 
transition between hover and wing-borne flight. 

NOMENCLATURE 

A jet area 

BL butt line 

C D drag coefficient 

C L lift coefficient 

C M pitching-moment coefficient 

C T thrust coefficient 

C y Jet momentum coefficient 

c t section lift coefficient 

Cp pressure coefficient 

c* section jet momentum coefficient 

D plate or Jet diameter 

d Jet diameter 

EVD elementary vortex distribution 

FS fuselage station 

H,h height 

hy upflow jet thickness 

I turbulent intensity 
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K h factor for effect of forward speed on 

suckdown 
L lift 

l stagnation line 

LID lift improvement device 

M Mach number 

MPR nozzle pressure ratio 

P,p pressure 

q dynamic pressure 

R radius, radial distance 

R n Reynolds number 

S wing area 

S’ distance between jet impingement points 

on a ground plane 
STOL short takeoff and landing 

STOVL short takeoff and vertical landing 

T thrust, temperature 

V velocity 

V e effective velocity ratio, /q/qj 

VSTOL vertical or short takeoff and landing 
WL water line 

X,Y,Z Cartesian coordinates 
a angle-of-attack 

8 sideslip angle 

AL aero/propulsion lift-interference 

increment 

AM aero/propulsion lift-interference 

increment 

AT inlet temperature rise, T^ - T Q 

5 deflection angle 

5’ section deflection angle 

<t> velocity potential 

$' computational polar angle measured in 

the ground plane about a jet impinge- 
ment point referenced to the line 
joining two jet impingement points 
P , o polar coordinates 

Subscripts 

e effective or exit 

edg edge 

ej jet exit 

F fountain 

f flap or fountain 

h height or hover 

i inlet or point of maximum rate of change 

of decay parameter 
in inlet face 

J jet 

L LID 

max maximum 

n nozzle 

N normal 

o ambient 

pi, p2 wall jet 

t total 

und under surface 


up up from ground plane to undersurface 

V ground vortex 

w wake 

x x direction 

® free stream 

Superscript 

’ section 

MULTI -JET OUT-OF -GROUND EFFECT HOVER 
INTERACTIONS 


During hover, there is a base Joss due to 
interactions between the lifting jets and the 
lower surface of the aircraft which results in a 
distribution of induced suction pressures which 
produce a lift loss. Additional performance 
losses include inlet-flow distortion, hot-gas 
ingestion, hot-day conditions, control bleed, 
internal nozzle flow, thrust vectoring, and 
static ground effect. There are many items 
related to the details of the aircraft design 
which determine the magnitude of the losses. 

Even though the sum of these losses may be only 
a few percent of rated thrust, an accurate 
knowledge of each is required to make a realis- 
tic estimate of the aircraft’s performance. An 
error of as little as 3% in thrust would reduce 
the gross weight, which in turn would reduce the 
fuel capacity, and hence the design range, by as 
much as 10$. 

The classic form of an axisymmetric free 
jet with a uniform nozzle exit dynamic pressure 
profile and low turbulence is sketched in 
Fig. 1. In the jet there are two regions of 
flow: a) the short potential-core region (up to 

six nozzle diameters long) which has a conical 
shape and a uniform velocity profile, and b) the 
fully turbulent region. In an early program 
intended to evaluate base losses during hover, a 
NASA investigation (17) evaluated the effect of 
the character of the jet on induced lift loss. 
The results showed that a relation exists 
between the lift loss and the rate of decay of 
nozzle dynamic pressure. 

Both of these parameters are functions of 
the amount of air entrained into the jet and the 
proximity of the entrainment to the plate. A 
correlation between the slope of the lift-loss 
curve and a parameter indicative of the dynamic- 
pressure decay was developed. As indicated in 
Fig. 2, this parameter is the maximum slope of 
dynamic-pressure decay divided by the distance 
downstream where that slope occurs. The data of 
(17) were used to obtain the correlation pre- 
sented in Fig. 2. 
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where 

S = planform area 
A = total Jet-exit area 

d 0 = diameter of an equivalent single nozzle 
having an area equal to the total of the 
areas of the multiple nozzles 

3qx/(P n - p)\ 

1 = maximum rate of decay of 

6 max dynamic pressure 


(x/dg)! = downstream distance at which the 

dynamic pressure decay rate is maximum 

It should be pointed out, however, that the 
dynamic-decay parameter is quite sensitive to 
the details of how the data are faired. Since 
the magnitude of the lift loss is usually a few 
percent of thrust, a potential estimation error 
of 10t amounts to only a few tenths of a percent 
of the thrust and is not very serious. The 
solid data symbols in Fig. 2 were obtained from 
a large-scale test which used a J-85 turbojet 
engine to provide a hot Jet (approximately 
1000° F, and pressure ratios up to 1.73). The 
resultant lift-loss data agree with the 
correlation derived from small-scale data. 

Subsequent investigations by Lockheed- 
Georgia (18), under contract to NASA on configu- 
rations derived from the XV-4B airplane, show 
that this type of correlation can be suitable 
for additional conf igurations. However, it was 
determined that the effect of jet-pressure ratio 
needed to be included in Eq. 1 by multiplying 
the right side of the equation by (pn/p)“^*^, 
and by changing the constant from -0.009 
to -0.016. 

A more direct, easier-to-use method for 
estimating these hover lift losses was developed 
by McDonnell Aircraft (see Equation in Fig. 33 
from (2)). Correlation of data from various 
single- and multiple-jet configurations resulted 
in the following expression 


where P pt is the total perimeter of all of the 
jets in the configuration. It is noted that 
Eq. 2 implicitly accounts for the higher decay 
rate of multiple-jet conf igurat ions in terms of 
equivalent Jet diameter, but does not account 
for higher decay rates caused by jet-exit condi- 
tions involving high entrainment rates. If 
higher than normal turbulence levels and decay 
rates are involved, Eq. 1 should be used. 

MULTI-JET IN-GROUND EFFECT HOVER INTERACTIONS 


During hover close to the ground, the 
lifting Jets impinge on the ground and form wall 
jets flowing radially outward from their 
impingement points. These outward-flowing wall 
jets entrain air beneath the configuration and 
lower the pressure on the bottom surface, caus- 
ing a jet-induced down load or suckdown . When 
the wall Jets from adjacent jets meet, an upwash 
or fountain flow is created between them. If 
there are only two jets, this upwash takes the 
form of a fan-shaped fountain as shown in 
Fig. 3. If there are more than two jets, a fan- 
shaped fountain is formed between each adjacent 
pair of jets, and a "fountain core" is formed 
near the centroid of the area enclosed by the 
jets where the fountain fans from each pair of 
Jets meet. When these fountain flows impinge on 
the lower surface of the configuration, their 
momentum imparts an upward force, or lift, that 
tends to offset the suckdown created by the 
entrainment action of the wall jets. 

The flows involved are highly turbulent and 
are not amenable to potential-flow analysis. 
Consequently , the only methods available for 
estimating their effects on the configuration 
(7,19) are empirically based. Both methods use 
essentially the same data base for multiple-jet 
configurations as does Wyatt's method (20) for 
estimating the suckdown and Yen’s method (21), 
or a similar approach, for estimating the foun- 
tain effects. Wyatt (20) developed an empirical 
method for estimating the suckdown for single- 
jet configurations that has become widely 
accepted. References (7) and (19) use slight 
modifications of Wyatt's method to estimate the 
suckdown for multiple- jet configurations; refer- 
ence (19) calculates the suckdown of the indi- 
vidual jets and sums them, and (7) uses Wyatt to 
calculate the suckdown due to an equivalent 
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single jet having the total area and thrust of 
all the Jets of the configuration. 

Yen’s method (21) for estimating the lift 
due to the impingement of the fountain basically 
calculates the flow in the sector of the wall 
Jet from each Jet that impinges on the lower 
surface of the configuration. The momentum of 
this flow, when it reaches the lower surface of 
the configuration, is estimated by applying 
empirically derived factors to account for the 
velocity decay in the flow along the ground 
(wall jet) and in the fountain. Unfortunately, 
the sum of the suckdown based on Wyatt’s method, 
and the fountain lift based on Yen's approach, 
do not equal the net induced lift for most 
multijet configurations. In fact, in some 
cases, as shown in Fig. 4, the twin- jet configu- 
ration has more net lift loss than the single- 
jet configuration. Figure 4 is from (22), which 
presents data on a twin-jet configuration, and 
on a single-jet configuration using one of the 
jets and a planform of half the area of the 
twin- jet configuration (thus maintaining the 
same planform-to-jet area ratio). The results 
indicate that, in addition to the fountain 
effect, an additional suckdown can be generated 
on multijet configurations. 

In (7) the increment of additional suckdown 
experienced by multiple jet configurations was 
extracted from the data base by subtracting the 
equivalent single- Jet suckdown (based on Wyatt) 
and the fountain increment (based on Yen) from 
the data. Expressions for estimating this 
multiple jet suckdown were developed from an 
empirical correlation of this multiple Jet 
suckdown with the geometric characteristics of 
the configuration. As shown in Fig. 5, the 
multijet suckdown term can be larger than the 
fountain increment. The probable cause of this 
additional suckdown is shown in Fig. 6 (from 
(23)). Vortex-like flows are formed between the 
fountain and each of the adjacent jets. The 
pressure distributions show that, as expected, 
the impingement of the .fountain flow produces 
high lifting pressures on the center region of 
the plate between the jets, but the vortex-like 
flows between the fountain and the jets induce 
equally strong suction pressures. The estimated 
suckdown for a single-jet configuration with the 
same planform-to-jet-area ratio would correspond 
to a suction-pressure coefficient, averaged 
across the area of the plate, about equal to the 
outer contour line shown in Fig. 7 
(C p = -0.004). Thus, both the lifting pressures 
and the additional suckdown pressures are much 
greater than the pressures induced on a single- 
jet configuration, and the question of whether 


there is a net lift gain or loss depends on 
which predominates. 

Unfortunately, at the time the methods of 
(7,19) were developed there were no other pres- 
sure data of the type shown in Fig. 7 that coula 
help as a guide in developing the method for 
estimating multiple-jet effects. More recently 
similar data (also for a twin- jet configuration) 
have been published in (24), which show the 
effects of height (Fig. 7) and jet deflection on 
these pressure distributions. Additional data 
of this type are needed to provide a better 
foundation for estimating the multi jet suckdown. 

The methods of (7,19) also include factors 
to account for the effects of lower-surface 
(fuselage) contour, and lift improvement devices 
(LIDs) to increase the fountain lift. These 
factors are based on a very limited data base, 
and the applicability of the methods to non- 
planar configurations is therefore limited. 
Nevertheless, the methods work well on configu- 
rations that fall within the data base. Fig- 
ure 8 shows good agreement between estimates by 
the method of (7), and data for a Harrier model 
(25), with and without LIDs. 

On the other hand, Fig. 9 shows a compari- 
son with more recently acquired data (26) for a 
configuration with a different jet-spacing-to- 
fuselage-width ratio. The method badly over- 
estimates the fountain contribution. There are 
other shortcomings in the available methods. 
Reference (26) presents a modification of 
Wyatt's method for single Jets to account for 
the effects of jet-pressure ratio that appeared 
to exist in much of the single- jet data base. 
However, more recent results from a full-scale 
investigation (27) of the suckdown induced on an 
8-foot-diameter plate by the exhaust of a J-97 
engine impinging on a ground plane, shows no 
effect of pressure ratio (Fig. 10). These 
results also show significantly more suckdown at 
low heights than was predicted by any of the 
available methods (Fig. 11). On the other hand, 
work by Lummus and Smith (22) has shown a small 
effect of both pressure ratio and jet turbulence 
(Fig. 12) on a twin-jet configuration. 

Before significant improvement can be made 
in the methods to estimate the suckdown for 
multiple- jet configurations, it will be neces- 
sary to understand and resolve these difficul- 
ties in predicting the suckdown for the simple 
single-jet case. The turbulence measured in 
(22) was the turbulence in the jet itself. 
However, it is the wall jet that is entraining 
air from between itself and the model, and it is 
probably the turbulence in the wall jet that is 
important. A significant factor in the turbu- 
lence of the jet leaving the model, and probably 


5 


in the turbulence of the wall Jet, is the ring 
vortices generated and traveling downstream in 
the shear layer of the jet (Fig. 13 — photo- 
graphed by Drubka and Nagib and taken from 
(28)). It is the growth of these vortices that 
entrains the surrounding air and causes the Jet 
to decay by transferring energy from the Jet to 
the entrained air, and the transition to turbu- 
lence and subsequent jet decay. 

Reference (29) has noted that when the jet 
impinges on the ground, these ring vortices 
expand outward with the wall jet. They probably 
play a major role in the entrainment of the air 
between the model and the ground. A better 
understanding is needed of the suckdown effects 
due to Jet- induced flows. These flows include 
the ring vortices, the turbulence of the flow in 
the Jet before it leaves the nozzle, and the 
wall-jet development and turbulence. A proper 
modeling of the major features of the flow, 
including the turbulence in both the free Jet 
and the wall Jet, will be required to develop an 
adequate method for predicting suckdown of 
single- jet configurations. 

Obtaining a better understanding of the 
flow and an adequate method for predicting the 
suckdown of single-jet configurations is a 
necessary first step to developing improved 
methods for predicting the ground effects on 
multiple jet configurations. Once these are in 
hand, the work of Kotansky (2), Saripalli (30), 
Jenkins and Hill (3D, Kind and Suthanthiran 
(32), Gilbert (33), Foley and Finley (34), Rizk 
and Menon (35), and Childs and Nixon (36) on 
fountain flows can be reexamined. 

The ring vortices created in the shear 
layer of the free jet (Fig. 13) that expand 
outward in the wall jet emanating from the 
impingement point may also be responsible for 
the extreme unsteadiness of the fountain flow. 
When these vortices meet at the base of the 
fountain (as depicted in Fig. 14), their veloc- 
ity components will not cancel unless they are 
exactly in phase. The relationship of these 
vortices to the fountain strength and unstead- 
iness should be examined. However, the primary 
need in connection with improving our ability to 
predict the ground effects on multi jet configu- 
rations is to obtain a better understanding of 
the multijet suckdown and the vortex-like flows 
between the fountain and the Jets that appar- 
ently generate it (Fig. 6). 

HOT GAS INGESTION 


When a jet VTOL configuration hovers in 
ground effect, the same fountain flows (Fig. 3) 


that can produce a favorable induced lift can 
also transport hot gases to the vicinity of an 
unfavorably placed inlet. The hot gas causes 
either elevated inlet temperatures or inlet 
temperature distortion across the engine fan 
face. Either of these factors will reduce the 
engine thrust. Alternate engine inlets may be 
needed to avoid the hot gas. The fountain flows 
have been the subject of numerous investigations 
((30-36) for example). The entrainment, veloc- 
ity decay, and spreading rates of the wall jet 
and fountain flows have been studied, and empir- 
ical methods for estimating these factors have 
been determined. Kotansky (2) summarized much 
of this work and presented the methods developed 
for estimating the location of the stagnation 
line on the ground between the jets (Fig. 15), 
the deflection angle of the upwash fountain if 
the jets are of unequal thrust, and the momentum 
of the flow in the fountain. 

Green and Zanine (37) have used the same 
data base to develop a method for estimating the 
inlet-temperature rise due to fountain flows. 
Their method uses the empirically determined 
velocity and temperature-decay rates developed 
in the same data base to trace the temperature 
decay from the impingement point, through the 
wall Jet and fountain, along the lower surface 
of the configuration to the vicinity of the 
inlet as depicted in Fig. 16. The method pre- 
dicts the stagnation lines on the ground and on 
the under surface of the configuration as shown 
in Fig. 17. The method adequately predicts the 
rise in inlet temperature as the aircraft 
approaches the ground (Fig. 18), and for some 
configurations predicts a reversal of tempera- 
ture rise at the lowest heights that appears to 
be confirmed by the data. 

The method of (37) does not include the 
effects of LIDs because the data base on these 
devices is too limited. It also applies only to 
hover in no-wind conditions. If there is a wind 
(or in STOL operation) the forward-flowing wall 
Jet will be turned back toward the aircraft in 
the form of the ground vortex, and this provides 
another path for hot gases to find their way to 
the vicinity of the inlet. 

The inlet temperature rise is often greater 
at forward speed than at hover. The inlet 
temperature rises with forward speed (Fig. 19) 
because, as the velocity increases, the leading 
edge of the ground-vortex flow field is moved 
closer to the configuration, and the flow path 
and time for mixing to reduce the temperature of 
the gas is shortened. Eventually, a maximum 
temperature is reached. Finally, as the veloc- 
ity is increased further, a point is reached 
where the hot gas cloud is blown below or behind 


6 



the inlet, and no inlet temperature rise is 
experienced. 

Rudimentary correlations of the maximum 
inlet temperature rise experienced are presented 
in (38,39). For configurations arranged to 
minimize ingestion, the height of the inlet 
above the ground was found to be the primary 
parameter. It was found that the maximum tem- 
perature rise tended to be inversely propor- 
tional to the square of the inlet height-to- 
diameter ratio. These results are based on a 
very small data base and should be used with 
caution. Most of the data did not go to high 
enough velocities to establish a correlation of 
the velocity needed to avoid ingestion, but a 
method of estimating this speed based on the 
effect of forward velocity on the depth and 
forward projection of the ground vortex flow 
field is presented in (39). Improvements in the 
ability to estimate inlet temperature rise will 
depend on developing better methods to predict 
both the fountain flow field and the ground- 
vortex flow field. 

JET FLAP IN- AND OUT-OF-GROUND EFFECT 


Propulsive lift began in the early days of 
aviation when the propeller slipstream flowing 
over the wing was used to attain increased 
lift. Configuration concepts developed include 
propeller slipstream, internally blown Jet flap, 
upper surface blown flaps, and externally blown 
flaps. Along with the concept development and 
accompanying experimental work (40,41), a number 
of theoretical and computational methods have 
been developed to predict the aerodynamic char- 
acteristics of these configurations. 

The first theoretical method for a wing 
with powered lift was developed by Spence 
(42,43). This method gave the lift and pressure 
distribution on a two-dimensional wing with a 
jet flap. Spence used a vortex sheet whose 
vorticity depended on the section momentum 
coefficient and curvature of the sheet to repre- 
sent the thin jet wake behind the wing. He 
obtained the general results shown in Fig. 20. 
The lift results in Spence's method are in the 
form of two equations — one for lift-curve slope 
and the other for lift due to jet deflection. 
Both equations contain terms that depend on the 
section-momentum coefficient. The lift-curve 
slope, for example, is (2w) modified by two 
additional terms that account for section- 
momentum coefficient effects. The pressure 
difference between the upper and lower surface 
is a function of jet-deflection angle, section- 
momentum coefficient, and chordwise location. 


In the pressure coefficient plot in Fig. 20, the 
computed pressure difference has been combined 
with the pressure distributions over an ellipse 
to obtain the total pressure distribution for an 
elliptical airfoil with a jet flap. This formu- 
lation is valid up to a point near the trailing 
edge as shown in Fig. 20. 

Based on Spence's vortex-sheet approxima- 
tion for the two-dimensional jet wake, a number 
of three-dimensional methods for jet flaps were 
then developed. The earliest three-dimensional 
Jet-flap methods, which included those of 
Maskell and Spence (44), Hartunian (45), and Das 
(46), assumed a large aspect ratio and an ellip- 
tical spanload distribution. More general 
three-dimensional jet-flap methods have been 
developed from methods developed for conven- 
tional wings (Fig. 21) — that is, they assumed 
potential flow and used numerical integration of 
influence coefficient matrices to obtain a 
solution. The methods extended conventional 
theory by including Spence's relationship for 
the high-momentum jet wake. This implies the 
assumptions of thin jet sheet and modest wake- 
deflection angles. Some of these methods 
include lifting-line methods by Lissaman (47) 
and by Lopez and Shen (48) and lifting-surface 
methods by Goldhammer, Lopez, and Shen (49). In 
addition, Hackett and Lyman (50) developed an 
equivalent mechanical-flap method. 

Reference (41) presents many comparisons 
between these theories and experimental data. 

One example is presented in Fig. 22 for a model 
with a jet-blown flap over the inboard two- 
thirds of the span, a 45° flap deflection, and 
jet-momentum coefficient values ranging from 0.5 
to 4.0. The lift, drag, and pitching-moment 
coefficient results were obtained experimentally 
using this model, and the corresponding calcu- 
lated values were obtained using the elementary 
vortex distribution (EVD) lifting-surface theory 
(49). The data show generally very good agree- 
ment between theory and experiment, especially 
for the lift and pitching-moment coefficients. 

It should be noted that the profile drag has not 
been included in the EVD calculations; adding it 
would bring the calculated drag into closer 
agreement with the experimental values. 

An example of an application of the jet- 
flap theory to a configuration for which the 
assumption of a thin jet exhausting near the 
trailing edge of the wing was violated is pre- 
sented in Fig. 23. The experimental data, pre- 
sented in (51), are for an external-blown flap 
configuration having four engines whose exhausts 
impinge on a double-slotted flap. Calculations 
were made for this configuration using experi- 
mentally measured distributions of the exhaust 
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deflection angle and momentum coefficient at the 
flap trailing edge (shown on the right in 
Fig. 23). Such calculations might not be prac- 
tical in the general situation, of course, since 
this information would not ordinarily be avail- 
able. However, this example shows that jet-flap 
methods can give useful results even though the 
powered-lift system and resulting flow pattern 
were more complicated than a simple Jet flap. 

The data show that the maximum and minimum 
deflection angles occurred on approximately the 
centerlines of the inboard and outboard engines, 
respectively, whereas the peak momentum coeffi- 
cient occurred between the centerlines. 

The results on the left in Fig. 23 show 
that the two jet-flap theories — Lissaman (47) 
and EVD (49) — give fairly good approximations of 
lift coefficient for momentum-coefficient values 
of 0 and 2, but not at a momentum coefficient of 
4. The EVD program underpredicted the lift, and 
the Lissaman program overpredicted it. At least 
part of this discrepancy could be attributed to 
the fact that the flap, rather than having a 
thin jet sheet emanating from its trailing edge, 
was operating within a relatively thick jet. As 
a result, noticeable Jet interference effects 
come into play which are not accounted for by 
the jet-flap theory. Further improvements are 
needed to deal with complete configurations 
including all of their complexities. Ultimate 
solutions for propulsive-lift configurations may 
require solution the use of Navier-Stokes 
equations. 

JET STOVL CONFIGURATIONS IN STOL OPERATION 


The prediction of the aerodynamic forces 
and moments experienced by jet- and fan-powered 
STOVL aircraft operating in ground proximity at 
forward speed (STOL operation) requires modeling 
a complex set of jet-induced flows in order to 
cover the entire height and speed range. These 
induced flows include the flow fields causing 
the induced lift losses in hover (in- and out- 
of-ground effect), the jet/free-stream interac- 
tion at forward speed, and the ground vortex 
created by the interaction of the free stream 
with the wall jet flowing forward from the Jet 
impingement point. Any method for predicting 
the aerodynamics of STOVL aircraft in STOL 
operation must reduce to a method which accounts 
only for the Jet- induced effects at forward 
speeds. 

An empirical method for predicting the 
effects of ground proximity in STOL operation 
was developed in (26), and is also available in 
(6). Reference (26) modifies the method of (7) 


for hover ground effects to account for the 
effects of forward speed, and modifies the 
method of section 2.2.1 of (4) (also in (54)) 
for the jet-induced effects in transition out- 
of-ground effects to account for the effects of 
ground proximity. To these is added an empiri- 
cal method, developed in (26), for estimating 
the ground-vortex effects, and a fourth term to 
account for forward-speed effects on the foun- 
tain contribution. The terms involved are shown 
schematically in Fig. 24. 

A key factor in the ground effects experi- 
enced in STOL operation is the location and 
strength of the ground vortex (Fig. 25). This 
ground vortex is formed where the forward flow- 
ing wall Jet is opposed by the free-stream flow 
and turned back on itself. This ground vortex 
induces suction pressures on the ground (and on 
the lower surface of the aircraft) as shown at 
the top of Fig. 26. An empirical correlation of 
the location of the ground vortex was developed 
in (26) (see also, Fig. 26), and used in the 
method for estimating the lift loss induced by 
the ground vortex, and for estimating the effect 
of forward speed on the hover suckdown term. 

As discussed above, and shown in Fig. 24, 
the method of (26) for estimating the lift and 
moment in-ground effect is made up of four 
terms. The ground-vortex term (L/T) v and the 
hover suckdown term (L/T) h induce a negative 
lift; the jet-wake term (L/T) w also usually 
induces a lift loss out-of-ground effect, but 
the proximity of the ground reduces this loss; 
and the fountain term (L/T) f produces a favor- 
able or positive lift. 

There is a limited data base with which to 
evaluate the method, and most of it is contained 
or referenced in (26). Figures 27 (using exper- 
imental data from (26)) and 28 (using experimen- 
tal data from (53)) indicate how well the method 
reproduces some of the data on which it is 
based. The estimates of the effects of ground 
proximity in STOL operation will be most reli- 
able if hover and out-of-ground-effect, jet- 
induced increments are available on the configu- 
ration as a starting point. Additional discus- 
sion of pertinent methods may be found in (52). 

JET INDUCED EFFECTS IN TRANSITION 


The transition flight regime ranges between 
hover and wing-borne flight, and uses the combi- 
nation of jet lift and aerodynamic lift to 
support the aircraft weight. The lifting jets 
induce significant changes in the flow field in 
the vicinity of the aerodynamic lifting sur- 
faces. The interactions between the lift jets 
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and the Tree-stream flow field are due to jet 
blockage, entrainment, separation in the jet 
wake, and lift-induced vortices. These interac- 
tions generate a complex flow field which is 
extremely difficult to analyze. 

Existing prediction methods (2-8) usually 
consist of either empirical procedures or 
potential-flow methods with empirical correc- 
tions to account for many features of the lift- 
jet/free-stream interactions listed in the 
previous paragraph. The empirical procedures 
tend to be easier to use and are often used in 
the preliminary design stage to identify the 
configurations which are the most appropriate 
for further design effort. The latter group of 
methods tend to require input data which are 
very time-consuming to prepare, and can require 
large amounts of computer time. As a result, 
they may be used for more detailed configuration 
evaluation; however, at their present state of 
development they are not consistently reliable 
nor are they adequately validated, and they 
cannot be counted on to provide more accurate 
results than the empirical methods. 

The most commonly used empirical method for 
estimating the jet- induced effects on VSTOL 
configurations in the transition flight regime 
is the one developed by Kuhn (in section 2.2.2 
of (4), and also in (54)). A less general 
version of this method which is suitable for a 
simplified preliminary design was presented as a 
computer program with a listing and sample input 
and output data (5). In either version of this 
method, its application is limited to configura- 
tions which are consistent with the data base 
used to develop the method. The approach used 
in this preliminary design method is outlined in 
the following paragraph. 

First, the pressure coefficient data 
obtained by Fearn (55,56) in an experimental 
investigation were used to develop the empirical 
formulation. These data were obtained for 
effective velocity ratios (V 0 ) ranging from 0.10 
to 0.45 in an experiment which had a 4-inch- 
diameter jet exhausting normally to the free 
stream, and from a large flat plate. The data 
were presented in pressure coefficients which 
represent the jet-induced increment. Equations 
were fit to represent the measured pressure 
distributions as a function of the effective 
velocity ratio, as well as the longitudinal and 
lateral distances from the jet exit. Next, 
these relations were integrated over the plan- 
form of the configuration of interest relative 
to the nozzle locations to estimate the jet- 
induced aerodynamics in transition flight. The 
results were then compared with available exper- 
imental data for configurations which ranged 


from simple flat plates (Figs. 29 and 30) to 
more realistic aircraft configurations 
(Fig. 31). These comparisons provided a basis 
for refining the estimates to account for con- 
figuration variables such as jet-nozzle locatic 
relative to the wing and fuselage, planform- 
to-jet exit area, jet deflection angle, and 
fuselage lower surface contour. These compari- 
sons also identified limitations of the met 'od 
(Fig. 32). Configurations with a jet located 
near the wing or other lifting-surface trailing 
edge have additional jet-flap-like induced lift 
which is inadequately accounted for by this 
method. These effects were discussed earlier in 
the section entitled "Jet Flap In- and Out-of- 
Ground Effect," and may be accounted for by any 
of the several extensions of Spence's jet-flap 
theory (40-47). 

Several potential-flow, vortex-lattice or 
panel-method-based schemes have beer developed 
which account for Jet-induced interactions. A 
comparison of five production-surface panel 
methods (57) demonstrated that good agreement 
for conventional aircraft aerodynamics can be 
achieved between panel-method results and 
experimental data. One of the earliest compre- 
hensive methods which used a vortex-lattice 
aircraft model (58) and a source-doublet Jet 
model was developed by Wooler (1,59). This 
method provided good comparisons between calcu- 
lated and experimental results for a limited 
number of VSTOL aircraft configurations. 

Examples of other methods have shown useful 
results (60-64). 

The complexity of predicting propulsive- 
induced effects on VSTOL aircraft prompted 
Wooler to use a modular approach (Fig. 33). The 
jet model is incompressible and neglects viscous 
effects other than the entrainment caused by 
potential flow sinks or doublets. The entrain- 
ment of free-stream flow into the jet and the 
pressure force on the jet boundary govern the 
equations of motion of the jet. The entrainment 
parameters are obtained from experimental 
data. The variation of the jet cross section is 
established, based on experimental observations, 
by a circular Jet transforming into an 
ellipse. The induced-velocity flow field due to 
the jet interference is obtained from two 
singularities: a) a uniform sink distribution 

on axes normal to the free stream at discrete 
locations along the jet centerline represent the 
entrained flow, and b) a doublet distribution 
along the jet centerline represents the blockage 
and Jet-induced circulation effects. Additional 
procedures are developed to represent 
nonaxisymmetric Jets, jet pairs, and a jet in a 
nonuniform free stream. 
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Another module evaluates the jet-induced 
forces and moments using a version of the 
vortex-lattice method from (58) to represent 
lifting planforms. The method also gives good 
agreement with experimental data for wing-body 
combinations by including the planform of the 
body in addition to the wing. Power-induced 
aerodynamic characteristics are evaluated by 
using a propulsion-induced camber distribution 
on the planform to satisfy the flow-tangency 
condition at the three-quarter chord of the 
horseshoe vortices of the lattice representa- 
tion. This wing-body program was used to 
improve the method of (1), which represented 
only the wing. 

Experimental data (65,66) from the model 
shown in Fig. 34 were used to evaluate these 
methods. These data were obtained in conjunc- 
tion with the research program (1). Using the 
method of (1), Mineck (67) obtained good agree- 
ment with experimental data from the configura- 
tion with the nozzles in the aft position 
(Fig. 35). For this case the jet-induced 
increase in wing lift (i.e., jet-flap effect) 
dominates. In contrast, poor agreement was 
obtained with the nozzles in the forward 
position (Fig. 35). The improved Wooler method 
of (59) was also used to evaluate these cases 
(Fig. 36). The wing-alone calculation based on 
the original method provided poor agreement for 
the forward jet position. To represent the 
presence of the large nacelles and the body, a 
vortex-lattice method (58) was used to include 
them in the planform for the improved Wooler 
method. The improved correlation is shown in 
Fig. 36. These results indicate that this 
improved method (59) by Wooler is suitable for 
propulsion-induced aerodynamic forces and 
moments at the conceptual and preliminary design 
stage. The output from the Jet- induced flow- 
field module can be used with a panel method for 
more detailed analysis. 

In an early application of panel methods to 
model V/STOL aero/propulsion interactions, 
Rubbert (68) used the Boeing TEA-230 panel 
method to represent a jet in a crossflow. The 
specification of the Neumann jet-boundary condi- 
tions proved to be difficult because the 
detailed distribution of jet entrainment around 
the curved, three dimensional lifting jet was 
unknown. A later attempt by Siclari (64) at 
Grumman was more successful. This method used 
empirical and theoretical data that describe the 
development of free jets to specify the needed 
jet-boundary conditions. A set of mixing data 
was developed which simulated a) the viscous 
entrainment of ambient air for both free jets 
and wail jets; b) the free- jet shape in 


cross-flow to provide the jet-blocking effect; 
and c) the jet formation along the ground which 
interacts with other jets leading to the 
"fountain” effect. 

These Jet data were combined with an exist- 
ing wing-body panel method (Fig. 37). The usual 
zero-normal-velocity boundary was used for the 
aircraft; the normal inflow velocity, simulating 
mass entrainment, was used for the jet and 
ground plane. In ground effect, the fountain 
momentum and incidence are calculated sepa- 
rately, with empirical data used to calculate 
additional fountain forces and moments. One 
correlation with test data from the VAK-191 VTOL 
aircraft is presented in Fig. 38 for lift and 
pitching moment coefficients. These data are 
for a lift engine thrust deflection of 77.5° and 
a lift/cruise thrust deflection of 60° at an 
effective velocity ratio of approximately 
0.14. At an angle-of-attack of 0°, the computed 
power-off results agree well with the experimen- 
tal data, and the computed power-on lift- 
coefficient results agree well with experimental 
data. As in the computed results from Wooler' s 
original method by Mineck (67), the poor agree- 
ment for power-on pitching moment indicates a 
need to include a separated-flow representation 
in the wake portion of the jet. 

A further indication of the need for a 
method which accounts for the lifting- jet-wake- 
induced separation is shown in Fig. 39 from a 
V/STOL aircraft transition flight analysis 
(69). This analysis used the PANAIR panel 
method to represent the aircraft with Neumann 
boundary condition and a parabolized Navier- 
Stokes solution for the lifting jet with an 
entrainment boundary condition. The jet shape 
and entrainment were determined using the Adler- 
Baron Jet-in-a-crossflow method (70) and recent 
unpublished experimental data acquired by 
Adler. The method was applied to a simple high- 
wing aircraft configuration with a single lift- 
ing jet. Unfortunately, the calculated results 
(Fig. 39) using this analysis do not account for 
the separated wake due to the lifting jet. As a 
result, the propuls ion- induced lift loss is 
underestimated by about 20 % of the thrust. 

CONCLUDING REMARKS 


For over 30 years, VSTOL aircraft have been 
developed primarily from experimental investiga- 
tions because there are many unique aerodynamic/ 
propulsion-induced effects which occur with 
these vehicles. The existing data base has been 
used to develop some empirical and modified- 
linear, potential-flow methods for predicting 
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the induced aerodynamic effects. The adequacy 
of available methods is graphically portrayed in 
Fig. 40. For propulsion induced effects 
incurred out-of-ground effect, available methods 
are adequate in hover and moderately adequate 
for transition flight. For propuls ion- induced 
effects incurred in-ground effect the 
interactions are more complicated and available 
methods are inadequate. 

Improved prediction methods are needed to 
guide experimental programs and to reduce the 
testing requirements in the future, as well as 
to provide an understanding of the flow 
physics. The Navy V/STOL Aerodynamics, Stabil- 
ity and Control Manual (4,5) utilizes empirical 
methods, linear-inviscid methods, and linear- 
inviscid plus boundary-layer methods, and repre- 
sents the current state of the art in V/STOL 
computational methods. These methods are very 
useful for preliminary design of configurations 
which are consistent with the underlying 
empirical data base. These methods are usually 
not suitable for detailed aircraft design. 

Multiple jet aircraft out-of-ground effect 
hover interactions represent a lift loss of a 
few percent of the total thrust, and these 
interactions may be estimated adequately by 
empirical procedures to within a percent of the 
total thrust. In ground effect, the lift-jet/ 
ground interaction creates complex near-field 
fountain flows which are difficult to analyze 
for induced forces and moments and for hot-gas 
ingestion. Jet-flap theories provide useful 
preliminary design estimates, but they do not 
account for all of the major jet interference 
effects. For lifting Jet-induced effects in 
transition flight, the empirical methods of Kuhn 
(54) are limited to configurations which are 
consistent with the data base used to develop 
the methods. 

While existing empirical methods are 
inadequate for detailed aircraft design, results 
from current applications of Reynolds-averaged 
Navier-Stokes solutions (12-16,35,36) offer a 
promising approach. Initial parabolic Navier- 
Stokes applications (12,69) to STOVL aircraft 
configurations do not adequately handle Jet- 
wake-induced flow-separation effects. While 
initial elliptic Navier-Stokes solutions are 
demonstrating greatly improved results, there is 
a need for improved turbulence models for the 
propulsion-dominated flows. Based on current 
CFD activity, the Reynolds-averaged Navier- 
Stokes solutions should achieve major improve- 
ments in the next several years and lead to 
analysis methods which can estimate aero/ 
propulsion interactions in a manner suitable for 
detailed aircraft design applications. 
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TURBULENT 
MIXING BOUNDARY 



Fig. 1. Schematic of the decay and spread of the jet efflux with distance downstream from 
the nozzle exit. 



METHOD 1 



O SINGLE JET 
□ MULTIJET 
A MULTISLOT 


MODEL 
REF. 6 



(b) 



Fig. 3. Fountain flow developed between a pair 
of lifting Jets. 



Fig. 2. Correlation of induced loads with jet 
decay parameter (p h _/p = 2.0). 

^ t K 


Fig. 4. Some multiple- jet configurations can 
experience more suckdown than an equivalent 
single-jet configuration (22). 
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Fig. 5. Fountain lift and the additional multiple-jet suckdown generated on a simple twin- 
jet configuration (7). 



Fig. 6. Flow field and induced pressure distribution for a twin-jet configuration (23). 
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Fig. 7. Three-view configuration and geometry (a) and pressure contours (b) of a twin- jet 
configuration (24). 
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Fig. 8. Lift induced on a Harrier-type config- 
uration hovering in ground effect (25). 
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Fig. 10. Effect of nozzle pressure ratio on 
single- Jet suckdown (27). 



O EXPERIMENTAL DATA, (REF. 26) 



Fig. 9. Lift induced on a simple twin- Jet high- 
wing configuration hovering in ground effect 
( 6 ). 



HEIGHT ABOVE GROUND, H/d 

Fig. 11. Suckdown induced by full-scale jet 
engines (27) compared with Wyatt correction (20) 
based on small-scale data. 
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Fig. 12. Turbulence screen and nozzle pressure 
ratio effects on net induced lift for a twin- Jet 
configuration (27). 



Fig. 13. Shadowgraph of Jet (R n = 10,000) 
showing initial laminar jet, the formation of 
ring vortices in the shear layer, and the 
transition to turbulence (28). 



Fig. 14. Shear layer vortices may be the cause of fountain unsteadiness. 
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STAGNATION 



POINT POINT 


Fig. 15. A method for estimating the location and strength of the upwash fountain generated 
between impinging jets is presented in (2). 


INLET FACE 



Fig. 16. A method for estimating the inlet 
temperature rise due to fountain flows by 
tracing the flow from impingement to the inlet 
is presented in (37). 



Fig. 17. Stagnation lines on the ground and on 
the aircraft lower surface calculated by the 
method of (37). 
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FORWARD SPEED. V, knots 


Fig. 18. Inlet temperature rise for a Harrier- Fig. 19. Typical variation of inlet temperature 

type aircraft estimated by the method of (37). rise with forward speed. 



Fig. 20. Variation of lift and pressure coefficients according to Spence's two-dimensional 
Jet-flap theory (43). 
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• V2 0 - 0 INVISCIO, INCOMPRESSIBLE FLOW 

• LIFTING-LINE OR LIFTING-SURFACE INFLUENCE COEFFICIENT 
METHOOS 

• WING AND WAKE COUPLED BY SPENCE JET-FLAP EQUATIONS 



LIFTING LINE LIFTING SURFACE 

Fig. 21. Schematic of three-dimensional Jet-flap wing lifting line and lifting surface 
methods . 



Fig. 22. Comparisons of EVD Jet-flap theory with experimental data for two-thirds span jet- 
blown flap model. 
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EST. (REF. 26) 


Fig. 28. Comparison of estimated ground effects with experimental data for a Harrier-type 
configuration of (53). 


MANUAL PREDICTION, REF. 4 
EXPERIMENTAL DATA. REF. 5 


d- ■ 3.5 in. 


Fig. 29. Comparison of predicted and 
experimental data of clipped delta-wing 
configuration for a single-jet configuration, 
SJ/S = 0.024. 


Fig. 30. Comparison of predicted and 
experimental data of clipped delta-wing 
configuration, SJ/S = 0.024. a) Single-jet 
configuration, b) Four-jet configuration. 








v . 



Fig. 32. Comparison of predicted and experimental data of a twin-jet, high-wing 
configuration, a) Nozzles in the forward position, b) Nozzles in the aft position. 



• DOUBLETS 
> SINKS 

+ DOWNWASH CONTROL PTS 


Fig. 33. Wooler method for computing the Jet-induced effects on the wing using superposition 
of lifting-line wing and sink/doublet Jet representations. 
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Fig. 36. Comparison among Wooler's original 
method (1), Wooler's improved method (59), and 
experimental data (1, 65, and 66) obtained on 
the model in Fig. 34 with the nozzles in the 
Forward position. 


Fig. 37. Three-dimensional singularity-panel 
approach to VSTOL aircraft prediction methods 
(64). 




Fig. 38. Correlation between three-dimensional singularity-panel method and experimental 
data (64). 
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Fig. 40. Adequacy of available prediction methods; shaded area indicates percent of need 
covered; numbers indicate key references. 
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